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i ' Abstract 

We show that in Euclidean field theories that have bounce solutions, the 
bounce with the least action is the global minimum of the action in an open 
space of field configurations. A rigorous upper bound on the minimal bounce 
action can therefore be obtained by finite numerical methods. This sets a 
lower bound on the tunneling rate which, fortunately, is often the more inter- 
esting and useful bound. We introduce a notion of reduction which allows this 
bound to be computed with less effort by reducing complicated field theories 
to simpler ones. 
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1 Introduction 



Vacuum tunneling in field theories was treated semiclassically in ref. M. In the small h 
limit the tunneling rate per unit volume can be expanded as: 

£ = Aexp(=|S) x[l + 0(h)] . (1.1) 

The exponent is the Euclidean action of the saddle point configuration called the "bounce" . 
The symbol denotes the entire field configuration at the bounce. The prefactor A has 
the dimensions of (mass) d in a d dimensional theory and can be formally written as a 
ratio of two operator determinants. The tunneling rate is usually estimated by replacing 
A by m d , where m is a characteristic mass scale of the theory. The error in making this 
approximation is not significant compared to the error coming from the exponential term 
which is determined by numerically computing the bounce action S[4>] (a more accurate 
estimate may be obtained in some cases, see ref. 0). 

The determination of T/V is important in extensions of the standard model where 
the effective potential may have several inequivalent vacua and the realistic and viable 
vacuum may turn out to be only a local minimum of the potential. In that case, we would 
be living in a false vacuum now, and our very existence would imply some upper bound on 
the tunneling rate. If a lower bound on the tunneling rate can be obtained by analytical 
or numerical methods, then a comparison with the experimental upper bound may place 
strong constraints on the parameters of the theory. One is therefore often interested in 
finding an upper bound on the bounce action by analytical or numerical methods. 

Recently there has been a renewed interest in considering the tunneling efects in models 
of particle physics J|, [5|, ||, [7], §. The effect of such tunneling was considered in ref. ||] 
to constrain the parameter space of the existing models of dynamical supersymmetry 
breaking (DSB). Similar tunneling rate computations exist in the literature || [j], || in the 
context of the minimal supersymmetric standard model (MSSM), where charge and color 
breaking (CCB) vacuua exist ||. While analytical estimates of the bounce action can 
be made in special cases [[I], the realistic models of DSB and MSSM often fall in the 
intractable cases where one must resort to numerical methods. When several fields and 
coupling constants are involved, the numerical search for the bounce becomes long and 
time consuming. The search also needs to be repeated many times to map a significant 
region of the parameter space. New techniques for approximating the bounce action in a 
short time are therefore welcome. 
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In this paper we point out the limitations of existing techniques of approximation and 
indicate a direction for further simplifications and improvements. In particular we show 
that all finite numerical methods are useful only to find an upper bound on the bounce 
action. We then develop a notion of "reduction" that provides an intuitive and systematic 
approach to better techniques for obtaining this bound. 

2 The Bounce as a Local Minimum 

To begin with, let us consider a theory with m scalar fields in d dimensions. The Euclidean 
action is 

S = T + V; (2.2) 

where T = Ea=ilJd d x (§^) 2 and V = J d d xU (0i, fa..<j> m ) with i = l,2...d and a = 
1,2. ..m. Suppose the potential U((f>i,<f)2---4> m ) nas a local minimum | at <p a = and a 
global minimum at <p a = <P a - Then a bounce solution 4> a (x) must possess the following 
three properties: (i) it approaches the false vacuum a = at t — ► ±oo; (ii) it is a saddle 
point of the action with a single "direction of instability"; (iii) at the turning point (say 
at t — 0), all the generalized velocities are zero. 

It will be useful to write down the effect of the scale transformation (x — > Ax) on T 
and V. In d dimensions one has: T — > T\ = \ d ~ 2 T and V — > V\ = X d V. Since the bounce 
is an extremum of the action, we must have, ■ASx{4 > ) \\=i = 0, from which it follows that 

(d-2)r(0) + (d)y(0) = o. (2.3) 

From (0) one also finds that ^5 A (0)| A =i = (d-2)(d-3)T$) + d(d-l)V(0) < when 
d > 2. Therefore the direction of instability associated with the bounce has a component 
corresponding to the scale transformations for d > 2. 

In [[TIJ it was shown that in d > 2 dimensional theories with a single real scalar field 
(a = 1), the bounce solution with the least action is an 0(d) invariant field configuration 
with the boundary conditions: d ^ r ' | r=0 = and 0i(r)| r ^oo = 0, where r is the radial 
coordinate in R d . In this case, the task of finding the bounce action numerically is rather 
easy. One simply solves the Euclidian equations of motion from some initial point <f>i(0) 
and looks at the limiting value </>i(oo). For arbitrary values of 0i(O) the value of 0i(oo) 
is either an "overshoot" (0i(oo) > 0) or an "undershoot" (0i(oo) < 0). Since the correct 
2 By a suitable redefinition of the fields, the local minimum can be made to lie at the point 4> a = 0. 
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value of 4>i(0) (called the "escape point") must lie between two trial values which end in 
an overshoot and an undershoot, the search converges rapidly by bisections. 

The task is considerably more complicated in the case of theories with many scalar 
fields. The bracketing property of the overshoot and the undershoot, which is obviously 
valid only when the field is a real scalar quantity, is lost. Therefore one can not search 
for the escape point by integrating the equations of motion. Because the bounce is not a 
minimum of the action, but a saddle point, a simple minimization of the action would not 



work either. For instance, the method suggested in ref. flTT| of maximizing the action of a 
field configuration on the lattice with respect to scale transormations and then minimizing 
it with repect to random variations of the fields is unlikely to converge to a bounce solution 
because variations that are orthogonal to a scale transformation are hard to identify in an 
actual lattice computation. A remarkable simplification is however achieved by adopting 
the method of ref. which consists of adding to the Euclidean action other terms that 
have the following two properties: 

(i) they vanish at the bounce, 

(ii) they remove the instability associated with scaling, i.e, the "improved" action with the 
new terms is minimized with respect to all variations including the scale transformation. 

When terms like these are added to the action, a bounce appears as a local minimum 
(instead of a saddle point) of the improved action. An example of the improved action is 

S{<t> a ) = S(<f> a )+Z n a n \A\^ , (2.4) 

where A = (d — 2)T + (d)V and p n and a n are positive numbers for n = 1, 2... . For d > 2 



the new terms added to the action have the two properties mentioned above ||12|| . 

Although the method of improved action reduces the problem to a pure minimization, 
the bounce is only a local minimum in an infinite dimensional space. In the absence of 
further information, it is by no means obvious that minimizations by finite numerical 
methods actually yield a useful result. Typically, the discretization involved in the nu- 
merical methods generates spurious minima of the improved action. It is not clear if the 
result of such a search yields an action that is greater or less than the bounce action. 
This is a drawback since one is often interested in a rigorous upper bound on the bounce 
action. Although finite numerical methods will always have some limitations (which we 
discuss later) we will now show that they can be sufficient to obtain such a rigorous upper 
bound. The approach suggested by us will also lead to new avenues of simplification and 
refinement of numerical techniques. 
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3 The Spherical Bounce as a Global Minimum 



Let us define a spherical bounce as an 0(d) symmetric field configuration with properties 
(i) and (ii) of the bounce provided the variations are restricted to preserve the 0(d) 
invariance. Clearly the 0(d) invariance of the configuration implies that it has property 
(iii) of the bounce. Our approach rests on the following two observations about the 
spherical bounces. f\ 

Statement 1: A spherical bounce is also a true bounce. 

Statement 2: The spherical bounce with the least action is a global minimum of 
the action in an open space of field configurations which obey the boundary conditions 
appropriate for a bounce. 

The proof of Statement 2 is simpler for d > 2. We will discuss the cases d = 1, 2 later. 
Before we sketch the proofs, let us first point out some immediate gains that result from 
the knowledge of the space C\. At the heart of our aproach is the trivial point that it is 
always easier to bound a global minimum from above than a local minimum. 

There are two general limitations associated with finite numerical methods. Firstly, 
in multi-scalar theories there may be several bounce solutions. The tunneling rate is 
determined by the bounce (or bounces) with the least action. There is no numerical 
method that guarantees convergence to the bounce with the least action. Secondly, any 
method of discretizing the improved action (|2.4j ) introduces spurious local minima which 
are local minima of the discretized improved action with respect to a finite dimensional 
space of variations, but are not local minima when all continuous variations of the fields 
are considered. These spurious minima may convey no information about any bounce 
action. 

These limitations preclude the possibility of finding the least bounce action to within 
controllable numerical errors in multi-scalar theories. However, as we have emphasized 
earlier, in many cases it is useful to find just an upper bound on the least bounce action. 
The knowledge of the space C\ allows one to approximate the bounce by a field configu- 
ration that lies exactly on the space C\ and provides a rigorous upper bound on the least 
spherical bounce action (which itself bounds from above the least bounce action). As it 
happens the space C\ is quite simple and it is easy to find field configurations lying on it. 
We now present the proofs of Statements 1 and 2. The latter also defines the space C\. 
3 This approach is closely related to the one suggested by A. Wipf 0. 
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Proof of Statement 1: This statement is in fact a consequence of the principle of 
symmetric criticality ||13|| . Writing the action of (|2.3| ) in polar coordinates we have 



where are the angular coordinates and are the measures corresponding to the angular 
gradients The quantity S 2 is positive definite and is minimized on the space of 0(d) 
invariant configurations. Therefore a spherical bounce is an extremum with respect to 
the 0(d) breaking variations too and no new negative eigenvalue is added to the Hessian 
by considering the 0(d) breaking variations. That is, the spherical bounce is also a 
true bounce. Q.E.D. 

Proof of Statement 2 (for d > 2): Consider the space of all 0(d) invariant field 
configurations that satisfy the correct boundary conditions and are extremized with re- 
spect to the scale transformations x — > Xx. Let us call this space C. All points in C obey 
( [2~3"1) . The space C naturally splits into two disconnected parts C — Co © Ci, where Co 
is the connected space containing the trivial solution <p a = where all fields assume their 
values at the the false vacuum at all times. At the trivial solution V = 0, and in a small 
ball around this solution V > because the false vacuum is a local minimum of energy. 
On the other hand from Q2.3|) it is clear that for points in C, V < 0. Thus Co consists of 
an isolated point. 

Points in C% have V < and T > 0. Every nontrivial configuration satisfying (|2.3|) 
belongs to C\. But the the point G which is the global minimum of the action in C\ satisfies 
the correct boundary conditions for the bounce, is a saddle point with a single direction 
of instability (the direction corresponding to the generator of scale transformations which 
is "orthogonal" to C\) and has = at t = by the 0(d) invariance. Therefore, by 
definition, it is the spherical bounce with the least action. Q.E.D. 

Comment 1: It is obvious that the action of an arbitrary field configuration in the 
space C is necessarily bounded below by either the action of the trivial solution Co or the 
action of a nontrivial solution of the equations of motion lying on C\. It is a remarkable 
property of the space C that it is always the latter. The essential point is this. In a 
discrete numerical method one allows for variations of the action over a finite dimensional 



S 



Si + S 2 , 




(3.5) 



5 



space only. If scale transformations are included in the variations and 0(d) invariance is 
imposed throughout the search, then one is assured that the "bounce" obtained by the 
numerical method lies on the space C\ and provides the required upper bound. 

Comment 2: The 0(d) symmetry is indispensible for any numerical method be- 
cause it reduces the number of variables drastically. However it is redundant from 
the point of view of defining the space C\. A milder symmetry like the Z 2 symmetry 
a (t,x) = 4>a{— t, x) is sufficient to ensure that G has property (iii) of the bounce. One 
can redefine the spaces C and C\ with the 0(d) symmetry replaced by the Z 2 symmetry. 
The corresponding point G is then truely the bounce with the minimum action. This is 
because, as one can easily show, the point G is a true bounce and, conversely, every true 
bounce has the Z 2 symmetry possessed by G. The latter assertion follows from observing 
that the equations of motion are invariant under t — > — t and imposing ^f-(t = 0) = 
implies that any solution (p a of the equations of motion also has this symmetry. 

Comment 3: There is more to gain from the above approach than the insight that 
numerical methods are useful only to find upper bounds on the bounce action. New 
and more effective techniques are suggested. This is motivated as follows. Underlying 
any numerical technique is a discretization of the action. For instance a straightforward 
discretization of the improved action with 0(d) symmetry is: 
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S(4>) = T(<f>) + V(cf>) + a\T(<f>)+2V(<f>)\ 
T(0) = OA^^n 3 ^ —ft J 

n=l a=l ZZA 

N 

V{<j>) = ^ d Y,n"U(^\<pt\..4t ] )- (3-6) 

71=1 

The radial coordinate in R d is discretized into N points and Q is the solid angle in d 
dimensions. The purpose of the discretization is to reduce the space of field variations to 
a finite dimensional space (to a space Cf here, which is homomorphic to R Nxm ). 

There is no particular reason to believe that the discretization in (|3.6j ) is the most 
convenient way of reducing the space of all variations to a finite dimensional space. Indeed, 
as we show below, tunneling in theories with many scalar fields can be extremely difficult to 
explore with the discretization of ( |3.6|) . Since the most significant part of the computations 
is to make the final point lie on the space C\ and not how the reduction to some space 
Cf is achieved, the natural question to ask is: why not implement the "reduction" of 
allowed variations at an early stage by reducing the complicated field theory to a simpler 
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one? The question itself suggests new and systematic ways of improving the numerical 
techniques. 



4 Methods of Reduction 

We will seek new methods of reduction. Our search is guided by the technical difficulties 
encountered with the discretization of ( |3.6| ). Let us briefly describe the essential features 
and the limitations of the reduction characterizing this discretization. We will call this 
method of reduction method 1. 



Method 1: With the discretization (|3.6|) one minimizes the action with the con- 
straint (d — 2)T + (d)V — 0. The constraint can be enforced by making the Lagrange 
multiplier a large. The result is a discrete trajectory 4> a (n) with n = 1,2...N, which 
can be extended to the continuous and piecewise linear function <fi a (r) given by </> a (V) = 
(1 — r + n)<j) a (n) + (r — n)<j) a (n + 1) for n < r < n + 1 and <p a (r) = <f) a (N) for r > N. 
When U is a polynomial in <p a , the functions T and V can be calculated to arbitrary 
accuracy using this function. The condition ( |2.3| ) may not be exactly satisfied by T and 
V obtained at this point. But one can always perform a scale transformation to satisfy 
( |2.3| ) exactly at the end. Thus the method of improved action is successful in bounding 
the bounce action provided one performs the necessary scale transformation at the end. 

In this case the search space for the minimization is iV x m dimensional, with m fixed 
by the theory. In an N x m dimensional Euclidean space, if one desires to reach within 
a distance of 5 from a minimum by a random iterative search (as suggested in ref. [0]), 
the number of steps required is about (^) Nm , where I is the maximum step length in 
any direction. One can improve this by using a multidimensional "greedy" minimization 
technique such as the Conjugate Gradient Method jnfl . The computational time typically 
grows as some power of N x m (depending on the complexity of the function). Also, the 
search space has many local minima with widely varying actions. Usually, one needs 
several iterations to arrive at a good one. Clearly the computation becomes harder as the 
number of fields m increases. Therefore a method of reduction where the computational 
complexity and time do not depend on m becomes desirable for theories with many fields. 
The alternative method of reduction presented below has precisely this virtue. 

Method 2: The action is given by ( |2.2|) . The false vauum is at <p a = 0. Choose a 
straight line, passing through the point <p a = 0, in the m dimensional Euclidean space 
R m with axes given by 0i,02--0m- The fields <p a can be constrained to take values only 
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on this straight line by putting <p a = y a (j) where y a are real numbers staisfying y£ — 1 
and is the reduced scalar field. Then the action is reduced to the reduced action S[4>} = 
J d d x 1/2(J^) 2 + U(<f)) which is a functional of a single real scalar field and the bounce 
action in this reduced theory is easily computed by the method of bisections described 
earlier. Note that the bounce exists if and only if the chosen straight line passes through 
some point (f) l a such that U (0) > U (</>^ . This method is an obvious reduction that 
intuitively seems to be a correct simplification. However, the relation of the "bounce" 
obtained in the reduced theory to the true bounce may not be immediately clear. But 
note that the action is extremized with respect to all variations that do not move the 
fields <f) a (x) out of the chosen straight line in field space. The crucial point is that scale 
transformations are included in these variations. Therefore the solution to the equation 



of motion satisfies ( |2.3| ) and by Statement 2 its action is rigorously an upper bound on 
the action of the bounce with least action in the full theory. The search time scales like 

N \ s ' when the search is made by the method of undershoots and overshoots, and 
has no dependence on m. The problem of multiple bounces and spurious minima does 
not usually arise. Also, the result is automatically a point on C\ without the need for 
adding new terms to the action as in (|2.4| ) . 

In practice one encounters polynomial potentials in <p a {x). It is often possible to find 
the position of the true vacuum or the saddle point in the potential between the true 
vacuum and the false vacuum. The lines joining these points to the false vacuum may 
yield fairly low values for the upper bound on the bounce action. One can do the search 
over several lines which can be judiciously chosen. 





















Method 







X 

Fig. 1. Comparison of Methods 1 and 2. The region above the dashed (solid) 
line is ruled out by method 1 (method 1). 

In Figure 1 we show a comparison of the performance of the above two methods in 
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a realistic DSB model [HJ whose vacuum stability has been studied in ref. ||. For our 
purpose it is sufficient to consider the scalar fields in the so called messenger sector of 
the theory which serves to communicate supersymmetry breaking to the fields in the 
standard model. For a natural range of parameters, the true vacuum in the messenger 
sector has been shown to be color breaking ||]. The relevant part of the action consists of 
the complex fields P, N, S, q and q. The last two fields carry color but the color index can 
be suppressed in this discussion. Apart from the usual kinetic energy terms, the action 
has the potential: 

V = ^-(\P\ 2 -\N\ 2 ) + (M 2 + Xj\S\ 2 )(\P\ 2 + \N\ 2 ) (4.7) 
+ K 2 \S\ 2 (\q\ 2 + \q\ 2 ) + \nqq + XS 2 + X X PN\ 2 + a 2 (\P\ 2 + \N\ 2 ) 2 , 

where g, A, Ai and k are positive coupling constants, Mf < and a 2 is of the order of 
10~ 2 . A false but phenomenologically viable vacuum exists at q = q = and \P\ = \N\ 
with S 7^ if Ai < (kXx)/(k + X) j|. We have computed the vacuum tunneling rate using 
both the methods described above. In Figure 1 we plot our results in the A, Ai plane for 
k = 1.0 and a 2 = 0.01. The region below the dotted line is allowed analytically, but the 
region above the dashed (solid) curve is ruled out by method 2 (method 1) because the 
lifetime of the false vacuum is shorter than the age of the universe. In this case a larger 
region in the parameter space is ruled out by using method 1, but the use of method 
2, which is much quicker to implement, can vastly reduce the region of parameter space 
that one must explore with method 1. There may also be theories where method 1 or its 
variations^ can provide superior bounds than method 2. A simple variation of method 2 
is obtained by replacing the straight line in R m by a curve. When several intuitive choices 
for the reduction are available, it is difficult to say which method will provide superior 
results. However method 2 is always the quickest one and should be used to rule out as 
much of the parameter space as possible. 



5 rf = 1,2 

We would like to remark on the cases d — 1 and d = 2. Statement 2 is in fact generically 
true for d = 2. To see this consider the deformation of the functional T in S to T{rf) 
defined by 

TM'ifM&r- (5 ' 8) 



We have suggested a different method of reduction in a related context in ref. (if 
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where rj is a positive real number. With this deformation the equation fl2.3| ) is modified 
for (d = 2) to 

ijT(rj) + 2V = 0. (5.9) 

Statement 2 is true as rj approaches zero from the positive direction. Thus G remains a 
saddle point in the sense that it is a maximum of the action with respect to some variation 
orthogonal to G\. However the negative eigenvalue of the Hessian |^ may approach zero 
as rj — > 0. But there is no underlying symmetry to make this situation generic. In other 
words, the eigenvalues of the Hessian are likely to be of the order of m 2 where m is some 
mass scale in the theory. Unless there is a symmetry to protect its smallness, an eigenvalue 
can not be made zero without fine tuning the parameters of the theory. 

The case d — 1 can be treated as follows. The equations of motion to be satisfied by 
the bounce are 

-w^w. (5 - 10) 

with the boundary conditions a (±oo) — > 0, ^f\ t =o = 0. Solutions to (|5.10| ) resemble 
the motion of a particle of unit mass in the potential —U(<f) a ). Let us define the surface 
D in the m dimensional space of the "fields" <p a as the surface where U = 0. The space 
D seperates as D = D Q © D\ with D given by the part that is connected to the point 
<f) a = 0. By definition D and Di are not connected in D. Consider the space of Z 2 
invariant trajectories <f> a (t) = 4> a (—t) with a (O) G D\, a (oo) — > D and satisfying the 
constraint U[(j) a (t)] > 0. The action of these trajectories is positive definite and has a 
global minimum, which, we claim, is generically a bounce. 

Let us briefly substantiate the claim which can also be proven along the lines of ref. 



17] . G is a minimum with respect to all allowed variations except variations that move 



the point <f) a (0) out of the surface D\. Generically, connected parts in D\ are m — 1 
dimensional surfaces. This leaves a single variation that moves a (O) out of D\. We 
need to show that there is a direction of instability associated with this variation. The 
entire trajectory is a solution of the equations of motion. By conservation of the "energy" 
E = J2™=i |(^") 2 ~ U(4>i,4>2---4> m ), the total energy is zero on this trajectory and the 
velocity is zero at the "highest point" (-jt(O) = 0). This implies that the generator of 
time translations has a node at t — 0. By time translational invariance of the bounce, 
the eigenvalue of the Hessian that corresponds to time translations is zero. Therefore 
there is a nodeless variation with a negative eigenvalue and G is a bounce. 
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